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Abstract 

ON . The /^-delayed a-decay of 16 N has been used to restrict the El fraction of the ground state 7- 

o 

transition in the astrophysically important 12 C(a,7) 16 reaction in several experiments including 



(N 



(N 



those performed at TRIUMF and several other laboratories. A review of the published mea- 
surements is given, and GEANT4 simulations and R-Matrix calculations are presented to further 
clarify the observed a spectra. A clear response-function, in the form of a low-energy tail from the 



^ ' scattering of a-particles in the catcher foil is observed in these simulations for any foil thickness. 

<D ' 

Contrary to claims in the literature, the simulations show that the TRIUMF measurement and 

O _ 

those performed at Yale and Mainz originate from the same underlying spectrum. The simulations 



suggests that the discrepancies between the Yale and TRIUMF final results can be attributed to 



^ • incorrect deconvolution methods applied in the former case. The simulations show in general that 

o: 

the form (width) of the spectrum is very sensitive to the catcher foil thickness. It is concluded that 

m 

the TRIUMF measurement most likely represents the currently closest approximation to the true 

Q ' /3-delayed a-decay spectrum of N. 
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I. INTRODUCTION 



The /3-delayed a-decay spectrum of 16 N has been used to restrict the El fraction of the 
ground state transition in 12 C(a,7) 16 0. This reaction has often been referred to as the "holy 
grail" of experimental Nuclear Astrophysics due both to the importance of the 12 C (a, 7) 16 O 
reaction in stellar evolution, and the experimental difficulty in measuring the cross section at 
low enough energies to be relevant for quasi static helium burning. Some possible assistance 
in determining the El ground state transition rate arises from the fact that, in the /3-delayed 
a-decay of 16 N, the subthreshold E x =7.117, J 7r =l~ state of 16 is populated with a yield of 
about 0.05 in the /3-decay of 16 N, while the unbound E x =9.60 MeV, J ?r =l~ state is populated 
with a yield of about 10 -5 . This leads to a relatively stronger influence of the subthreshold 
state on the /3-delayed a-decay spectrum of 16 N, particularly on its low energy part. As the 
ground state of 16 N has a spin and parity of J 7r =2 _ , only 16 states with J n =l~ and 3~ 
will contribute to allowed (3 a decay of 16 N, with no J 7r =3~ state being nearby the given 
energy range of the aforementioned a spectrum. Therefore the 3^=3^ contribution to the 
spectrum is expected to be small, but not necessarily entirely negligible, while the relevant 
J w =l~ part is dominant. 

While there is, of course, a true, underlying /5-delayed a-decay spectrum of 16 N, any 
measurement of the spectrum will lead to experiment-specific distortions on the spectrum, 
like those from detector resolution, collector foil effects and possible backgrounds. While 
the last of these can be removed by a correct subtraction, it is quite difficult to remove 
the former effects as they convolute the ideal spectrum with a usually energy- dependent 
response function. In this case each bin point of the measured spectrum will be mixed more 
or less with contributions from elsewhere in the spectrum. In general, deconvolutions 35[] 



are considered a difficult mathematical problem with uncertainties growing rapidly the more 
iteration steps are applied and typically they involve some judgment calls. We will show 
that minimizing the system response is the only practical solution and that all attempts of 
a deconvolution lead to failure. 

As response functions (likely) differ in technically different experiments, experimental 
spectra of the /3-delayed a-decay of 16 N may in principle look rather different in compar- 
ison without implying a major difference in the underlying spectrum. In a correct fitting 
procedure, a theoretical curve should be convoluted with an experimentally well known re- 
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sponse function, and then be compared with the data; the theoretical curve can then be 
varied to fit the data after convolution. The theoretical fits may then be compared among 
different experiments. Of course, the sensitivity and the response of an experiment should 
be estimated beforehand independent of the actual 16 N spectrum. 

There have been several measurements of the /3-delayed a-decay spectrum of 16 N 

These measurements will be discussed below and compared, with the 
measurement described in Ref. Q serving as a benchmark. In Sec. [TT] we will describe 
individual measurements and their respective methods to obtain their final results, in Sec. II III 
we will present GEANT IV simulations of the experiments, in Sec. II VI different measurements 
are compared to each other, in Sec. |V] the final conclusions will be drawn and in Sec. IVII 
some future possibilities regarding the /3-delayed a-decay spectrum of 16 N will be discussed. 
In particular, besides many other issues of more minor concern, it will be shown that 

• the spectrum obtained at Mainz University suffers from considerable distortion effects; 

• the method applied to deconvolute the spectra obtained at Yale University appears to 
be mathematically erroneous; 

• the measurement done in Argonne National Laboratory likely suffers from an unre- 
ported background. 

II. DISCUSSION OF MEASUREMENTS 
A. The Mainz spectrum 



ana 



In the late 1960's to early 1970's a group at Mainz University [36j] led by H. Waffler made 
a series of very high statistics measurements of the /3-delayed a spectrum of 16 N 
with the ultimate goal of detecting the parity-forbidden a decay of the Ea;=8.872 MeV J 7r =2~ 
state in 16 0. A gas cell containing 15 N gas was bombarded with a deuteron beam, producing 
16 N via the (d,p) reaction. The activated gas was pumped into (one or) two detection 
cells, where the a decay was registered by (four or) eight silicon detectors mounted on the 
walls of the (one or) two counting cells. The silicon detectors were separated from the 6-8 
Torr gas volume by 30 /ig/cm 2 collodion (C6H 7 N 2 . 5 Oio) foils. No coincidence requirements 
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G. 1: (color online) Comparison of the Mainz spectrum as calibrated in Ref. [6| (cross) and Ref. 
(box); the energy is displayed in the CM system. 

were placed on the events to prevent detection of a-particles possibly degraded by system 
response. 

In 1971 H. Waffler communicated about a quarter of the events from the then available 
spectrum 



lOj to both C.A Barnes (Caltech) and F. Barker (Australian National University) 
(about 3.2xl0 7 events). These letters contained a spectrum as counts/bin sorted by channel, 
and a very precise energy calibration of incompletely explained origin. The data sent were 
selected for the smallest /3-induced tail of low energy pulses, for which the silicon detector 
biases were lowered to 8 V to reduce the depletion depth in the silicon. It was confirmed, 
however, by H. Waffler and by the position of the a peak from the singly forbidden /3-decay 
to the E x =9.85 J 7r =2 + 16 state in this spectrum that no correction for the 30 ug/cm 2 

nn 

collodion foils has been included in the calibration [12J, [13J . The matter was also discussed 



in Ref. 



14[ | . This spectrum, though unpublished, has been widely used in comparison to 



later measurements of the /9-delayed a spectrum of 16 N (see, e.g. Ref. jSj) . In particular, 
6] a table of an energy calibrated spectrum of the Mainz data is published (Tab. 



in Ref. 

A.l) without an adjustment for the foil thickness 371] . Fig. [1] shows a comparison between 
the two calibrations. 
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B. Spectra from the Yale group 



1. France III et al. publications 

As there is a recent publication by the Yale/U. of Conn, group [6] about their second 
measurement of the /^-delayed a spectrum of 16 N, this work will be discussed here first among 
the two Yale measurements. However, the first measurement of this group j^], published 
in 1993, had basically the same features as the second experiment with the exception of a 
different data reduction from the raw (/^-efficiency corrected) spectrum to the final spectrum, 
and, as a result, a vastly different final spectrum (Sec. IIV A II) . 

In both measurements the d( 15 N,p) 16 N reaction was used for production with the 16 N 
escaping the deuterium target at high velocity. To catch the 16 N nuclei from the low energy 
part of the recoiling 16 N velocity distribution a tilted aluminum foil of 180 /ig/cm 2 thickness 
was used. After the collection period, the catcher foil was moved in front of a silicon counter 
array with /3-counters positioned behind the foil. A (3-a coincidence was used to discriminate 
against some kinds of background events. To obtain the spectrum that was deconvoluted 
later, a /3-efhciency correction was necessary, particularly for the higher energy a events 
where the coincident (3 energies are very low. Fig. [2] shows a comparison between the 
measured and the /3-efhciency-corrected a spectrum and the ratio of the two spectra. The 
correction rises with increasing energy to a factor of 4.5 and then falls for the highest energy 
point to about 2.5. Little is said in Refs. how this correction for the extremely low 

energy (3 particles at the high a-energy side was actually derived. With decreasing f3 energy 
one would expect, in general, a smoothly rising ratio, i.e. the inverted efficiency curve, and 
a possible detection threshold for the /3-rays. 



About 235,000 events were produced 



6| from the analysis of the time versus energy 



spectrum. The raw a spectrum (/^-efficiency corrected, 285,000 events) from this analysis 
(displayed in Ref. Q], Fig. 5.5) is shown in Figj3]in comparison with the published TRIUMF 
spectrum j^j]. It is quite obvious from Fig. [3] that the spectrum of Ref. j^, 6] is both shifted 
and broadened compared to the TRIUMF spectrum. Also the low energy plateau is about an 
order of magnitude higher than in the TRIUMF measurement. This is clearly an effect of the 



very thick collector foil compared with the one used at r 
principle, a careful simulation of the spectrum in Ref. 0, 



TilUMF (10 /ig/cm 2 of carbon). In 
6| might have produced information 
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FIG. 2: (color online) Upper Panel: Comparison between the measured spectrum in Refs. jf], 
f|(box) and the /^-corrected spectrum (cross). Lower Panel: Ratio of the two spectra. The energy 



is in the laboratory system. 



about the underlying /3-delayed a spectrum, though likely with a large uncertainty and 
considerable dependence on the model used, see Sec. IIII C 21 However, in Refs. (a, la] a 
different approach to the data reduction was chosen, described as follows: 



'The measured spectral line shape was corrected for distortions caused by 
the variability of our time and energy resolution, see Fig. 5. For a spectrum 
constant in energy the yield measured at each point in that spectrum is directly 
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FIG. 3: (color online) Comparison of the final TRIUMF |2] (bars) and the /^-efficiency corrected 
a-spectrum of Ref. [5[ (crosses). The spectrum of Ref. |5[ has been normalized to the height of 
the TRIUMF spectrum. The abscissa represents the ce-energy in the laboratory. 

proportional to the energy integration interval. This conclusion holds for a 
spectrum in any physical variable. In the case of this experiment, the data are 



integrated over time anc 
energy resolutions. [.. 



energy with the integration intervals being the time and 
The final spectrum, with these resolutions divided 



out, is shown in Fig. 6 [....].' 

Note that no reference to an experiment is given, in which a similar deconvolution method 
has been used. 

The word 'resolution' in regard to the energy is used here in a somewhat unusual way, 
as it clearly labels the energy-dependent energy loss through the aluminum foil. The same 
is true for the 'time resolution' which is just the time of flight between the collector foil 
and the a-particle detector, and not the experimental spread of this time-of-flight signal. 
The final spectrum shown in Refs. is indeed derived by division by this 'energy 

resolution' and this 'time resolution'. Questions of dimensionality, i.e. divisions by energy 
and time, and the obvious renormalization that must have been performed to conserve the 
total number of counts are not discussed in these references. 

The claim that an experimental convolution, i.e. the result of a bad resolution, can be 
undone by dividing by an arithmetic function 39] that is independent of the final spectrum is 
not correct. Such a general claim can be easily refuted by a counter-example. In many cases 
in experimental physics a <5-like function (e.g. a narrow line) is folded in first order due to 
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FIG. 4: (color online) A Monte Carlo simulation (1000000 events) of the /^-delayed a-decay spec- 
trum of 16 N as found in Ref. [2| (bar points) exposed to the energy loss as described in Refs. a, \6\ 
(cross points), and a comparison to the France III raw spectrum (asterisk points), appropriately 
matched in height. 



statistical processes in the event collection by a Gaussian convolution. No division by a finite 
number will recover the initial 5-like function from the Gaussian measured. The number of 
events will stay conserved by the convolution. If only a fraction of the Gaussian function is 
integrated, the number of counts will, of course, change depending on the integration interval. 
Such a change is neither linear, nor will the total number of counts ever be exceeded, however 
far one may choose the integration interval. A reasonable integration interval for such a peak 
will, of course, depend on the resolution, but is not identical to it. In fact, by this method 
of division, for a 5-shaped input function, a 5-shaped output function results as the input 
only needs to be multiplied by the constant resolution. I.e. no degradation or convolution 
of such a line will result. It may be argued that the situation is different for a continuous 
spectrum, where, indeed, the resolution may be energy dependent. However, that would 
make the response function dependent on the spectral form, which is not expected from a 
reasonable experiment. 

It is straightforward to simulate, how the bare energy loss from the Yale catcher foil will 
influence the /3-delayed a-decay spectrum of 16 N. The results of the Monte Carlo simulations 
are shown in Fig. |H It is obvious that the energy loss effect gives a good description of the 
high energy side of the France III spectrum but misses further response effects on the low 
energy side (see Sec. IIIII) . Therefore there is more to the response function than the simple 
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stopping power effect (see Sec IIII CT|) . 



The main contribution to the deconvolution by division in Ref. [6[, already shown to 
be incorrect, comes from the time of flight resolution as displayed in their Fig. 5. as it 
shows a steep energy dependence. This use of the time-of-flight spectrum has no physical 
justification, as (i) the time-of-flight measurement does not influence the pulse height in 
the energy spectrum; (ii) the time-of-flight through the foil is about 50 fs for a 1 MeV a- 
particle, therefore the implantation (start) position in the catcher foil has nothing to do with 
the observed time-of-flight. The increase in the time of flight for lower energy a-particles 
should follow 1 / \f~E~ which is clearly exceeded in Fig. 5 of Ref. [6|]. However, as the spectrum 



discussed in Refs. (4, 5, 6] is not derived from the time-of-flight information, but from the 
pulse height in the silicon detectors, the time-of-flight, not influenced by the catcher foil 
thickness, has nothing to do with the response to the foil, but only with some kinds of 
background discrimination applied. For example, events degraded in the silicon detectors 
will be removed by the time-of-flight information. 

Besides a correction by division in the yield, the position of the energy points (initial bins) 
has also been changed, i.e. shifted to higher energies. Of course, a correct deconvolution 
does cause a shift and yield a change of channels simultaneously, see Sec. IIII C 21 This is 
claimed to be done in the following way 

'The effective energy of the emerging a-particles for each data point was 



calculated using the expected variation 40] of the yield over the energy width of 



the catcher foil for each slice. [...] Note that due to fast variation in the yield, 
the effective a-particle energy is not the one due to a-particles emitted from the 
center of the catcher foil.' 

This remark appears to suggest that a preconceived knowledge of the final spectrum has 



been usee 



in Refs. 



to derive these corrections. It can now be checked, how far the energy shifts used 



lis 



6J are consistent with the spectra shown in Fig. HI To calculate the effective 



energy E e ff we solve numerically the integral 

fl A{E) E'w(E>)dE> 
£ff SL HE) AE>)dE> 

for the effective energy. Here E is the initial energy, A the foil thickness, and w(E) the 
/3-delayed a-distribution of 16 N using i?-matrix theory applied to the TRIUMF spectrum. 
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FIG. 5: (color online) Expected energy shift for the foil used in Refs. a, |6[ (line) using Eq. [TJ The 
original and the shifted spectrum are shown in Fig. ??. In comparison, the energy shift used in 
Refs. 0,0] is shown here (bars). The energy is in the laboratory system. 

For simplification, the minor change in stopping power over the integration integral is being 
ignored. To calculate the energy shift, we subtract the effective energy from the initial energy. 
Fig. [5] displays the energy shift found here and the one used in Refs. 

aa 

. Compared to 

the Monte Carlo calculations, taking also the energy dependence of the stopping power into 
account, the curve shown in Fig. slightly underestimates the energy shift encountered. 
While the absolute position of these curves on the abscissa may be disputed, it is clear that 



the energy dependence given 
linear relation used in Refs. 



3V the calculation shown differs greatly from the approximately 



, |6(]. The use of the France III final spectrum as input would 
change little, if employed in the simulations. In fact, it is the quickly changing yield over 
the energy loss range in the foil of the /3-delayed a spectrum that causes these quite distinct 
changes in the energy correction. It is apparent that the energy correction used in Refs. 
0, [f| is equivalent to taking the center of the catcher foil employing the nearly linear energy 
dependence of the stopping power despite the opposite claim in the above quotation. This 
leads on the low energy side of the spectrum to an underestimate of the necessary energy 



correction in Refs 



.a 



while the correction is considerably overdone on the high energy 



side. Most importantly, calculating 'the effective a-particle energy' requires a knowledge of 
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the true distribution beforehand. 



2. Discussion of Z.Zhao et al. 



As m en«oned above, in a _nent p re v 10US to Ref, U Q a Relayed o- S pee trUm 
was obtained at Yale University by largely the same group |3J]. The experimental method 
was approximately identical to the one described for France III, except that the a collection 
system was changed from upstream to downstream of the target. Therefore, a spectrum very 
similar to the one shown for France III was collected (See Fig. la of Ref. 3]). However, the 
deconvolution method was different. As this measurement has not been indicated as incorrect 
or superseded in the later publications by this group and has been used for comparison in a 
recent publication Q], it will be discussed in the subsequent section. 

In Ref. j^l an .R-matrix description of the spectrum measured is given, following the 
parametrization of Ref. [if]]. Beyond this, the authors present in their Fig. lb the "£=3 
component of our fit" and conclude that "only with the introduction of the small £=3 
component was it possible to reproduce the line shape of the unfolded spectrum". It has 
since been stated [l7] that the £=3 component of Ref. jsj] is not the result of a fit, but 
simply "ad-hoc" , i.e. that it is an invented assumption by the authors without any basis in 



theory 



4l|. 



To see, however, the results of a consistent analysis of the data of Ref. 
analyzed those data using a full i?-matrix fit, following the description of Ref. 



31 we have 
. For the 



best fit a value of x 2 /point=0.17 was found for the data of Ref. [3J. This unexpectedly low 
value is attributed mainly to the data points of the main peak of Fig. 1(b) of ^] which 
show smaller fluctuations than those of Fig. 1(a), a result not normally expected from the 
deconvolution of an experimental spectrum (see Sec. IIII C 2\i . To further understand these 
findings, the deconvolution process based on the thesis of Z. Zhao 15] will be subsequently 
discussed and the likely reason for the low xt presented. 

The first step of the data reduction, after a /^-efficiency correction (Sec. Ill B lj) . was to 
correct for the energy loss in the target thickness. For this, it appears that the maximum 
of the experimental spectrum was shifted to match the maximum of the Mainz spectrum 
without further corrections. 

Independent of the energy shift, the number of counts in each bin has also undergone a 
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data reduction procedure. As stated in Ref. [3|, the final result is expected to be similar 
to that of Ref. 1], which is referred to as the "zero target thickness" spectrum |42] - Since 
the spectrum measured was far from agreeing with the one of Ref. [l| , the authors chose to 
complete the unfolding of the data by finding a response function which is to be determined 
"semi-empirically" [3J . 

A response function, [G(E, E — E') x R(E)], was then constructed by using it to convolute 
the so-called "zero target thickness" spectrum to yield the experimental Yale spectrum 
exactly. The first factor, the function G(E, E — E') corresponds to the simple folding over 
the catcher thickness, examples of such a folding have been shown above (e.g. Sec. Ill B 11 
Fig. 0J. Because this folding is still inadequate, a further correction factor, R(E), the 
second factor above, was required. This multiplicative correction factor R(E) was then 
determined precisely by taking the ratio of the experimental curve to the curve convoluted 
with G(E, E — E'). Note that in this procedure, G(E, E — E') does not have to be realistic, as 
R(E) takes into account all shortcomings inherent in the process. In fact, G(E, E — E') can 
be set to 1.0 at all energies with the same final result. Again a simple arithmetic function is 
used to perform a deconvolution and the response function is not independent of the actual 
a-spectrum. 

The procedure described represents the 'semi-empirical' [3( determination of the response 
function. The circular logic of this data reduction is obvious: if the "folding" of the "zero 
thickness spectrum" , i.e. the spectrum of Ref. with G(E, E — E') x R(E) yields exactly 
the experimental spectrum of Ref. g then the "unfolding- of the spectrum of Ref. with 
G(E, E — E) x R(E) will yield exactly the spectrum of Ref. pj. Expressed differently: the 

1] within the energy region given in 



data of of Ref. 



reproduction of those of Ref. 
Ref. p|, albeit with larger error bars 43] and fewer total counts. This explains the major 
part of the low x 2 /point=0.17 in the total .R-matrix fit to these data. 

It may be argued that the identity derived above was intended and that the procedure was 
designed specifically to provide a means of extrapolating the correction factor R(E) down 
to the low energy region, and that data points above Eq,=1500 keV are thus irrelevant to 
the analysis. While such a method would remain unphysical, there is no evidence anywhere 
in Refs. (3), [ijj] to indicate such a possibility. Factually in any serious i?-matrix analysis, 
these points at higher energy are of utmost importance. As displayed in their Fig. lb, the 
.R-matrix analysis in Ref. [3( suggests indeed the use of the entire spectrum. 
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FIG. 6: The function R(E) as derived in Ref. [3j versus center-of-mass energies for the lower energy 
part of the spectrum. 

The extrapolation of R{E) to lower energies is not described in Ref. [3J]; in Ref. [15( the 
following quotation can be found (Sec. 4.1, p. 87): 

'Consequently, the unfolding procedure carries small uncertainty down to 1.5 
MeV and we further extrapolate it down to 1.1 MeV. Since R(E) reflects a tail 
from higher energies, the extrapolation is such that R{E) below 1.45 MeV is a 
constant value of 3.6.' 

With this statement the function R(E) for the low energy part of the spectrum is shown in 
Fig. [6j It is clear that this choice of an arbitrarily constant function R(E) after a steep in- 
crease in the "known" region is unjustified, even, if the deconvolution procedure was correct. 
More reasonably R{E) should follow some functional form which is essentially unknown be- 
low E=1.45 MeV; and could, in principle, even increase or decrease rapidly with increasing 
energy. However, R(E) is not a response function of the system by any conventional use of 



the word 44|, but an arithmetic divider relying on preconceived knowledge. Thus R(E) has 
no experimental basis. 

Applying the constant R{E) to the initial spectrum shows immediately that the low 
energy points are still above the final points, so that some different kind of 'deconvolution' 
step must have been involved. A suggestion of what was really done is given in the following 



quotation from Ref. 



(P-87): 
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'In fig. 43 we show the convolution of the theoretical curve corresponding 
to S^i=95 keV-barn with the response function, together with the experimental 
data.' 



Thus, the /3-delayed a-decay spectrum of 16 N as presented in Ref. [3| is, in the high 
energy region, a close reproduction of the data of Ref. while for the low energy points 
it closely follows a purely theoretical curve. It does not come as any surprise that the final 
result of Ref. [3j is indeed S,ei(300)= 95 keV b. 

It is also noteworthy that the fractional errors of the low energy data points of Fig. 
1(b) [3I are no larger than those of Fig. 1(a) in spite of the major unfolding procedure 
required. If Fig. 1(b) is converted to counts, some of the lowest energy points have about 20 
counts/channel for which statistical errors of at least ±25-30% are expected, even excluding 
errors from the unfolding procedure. However, e.g., for the energy point at 1.448 MeV, the 
yield is 20.2±2.5, i.e. with an error 0.56 of the expected square root of the number of counts. 
If realistic errors were taken into account, the errors on R(E) should be added in quadrature 
with those resulting from statistics, when the data reduction is done. 

There has been a suggestion that the errors in the /3-delayed a-decay spectrum of 16 N of 
Ref. jsj are correlated. It is, however, abundantly clear from the above procedure that the 
final errors should be uncorrelated, if the initial experimental errors are uncorrelated. The 
unfolding procedure is achieved simply by division of correction factors, a procedure that 
does not induce correlations between data points. There has, however, been a claim in Ref. 
[15I that the errors of the data points are correlated as the detector resolution is about 50 
keV, while the catcher thickness is up to 250 keV. Such a claim would be incorrect, unless 
some procedure had been employed which moves events among the channels after the data 
acquisition. There is no evidence of such data shuffling anywhere in Refs. |3|, ll5|] as each 
point in the spectrum is divided by an individual number labeled as the response correction. 
The low x 2 /point found in the .R-matrix fits is completely explained by the fact that the (3- 
delayed a-decay spectrum of 16 N of Ref. 3j is basically a reproduction of the high-statistics 
spectrum of Ref. [l| over most of the energy range, combined with a normalization to a 
theoretical curve for the lower energy data. 
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C. The Seattle spectrum 



In 1993/94 a measurement of the /3-delayed a-spectrum of 16 N was started at the Univer- 
sity of Washington in Seattle. While this measurement was never published by those who 
carried out the measurement, the data were, however, distributed for comparison and made 
available for publication by other groups. Particularly the spectrum is published in tabular 
form in Ref. [6j. The only, rather incomplete, information available about the experiment 
is from two annual research reports of 1994 [jj and 1995 [rs| . 

The experiment used a Ti 15 N target bombarded by a deuterium beam to produce 16 N 
which was partially captured in carbon foils of 10 or 20 fig/ cm 2 thickness. The activity was 
then transfered to a position between two silicon detectors, and a and 12 C particles were 
observed in coincidence, as in Ref. [sj. It is not clear whether the distributed spectrum was 
made with a 10 or 20 fig/ cm 2 thickness, or, if both spectra were combined. Other unknown 
information includes the method of energy calibration and the resolution of the detectors. 
It is, however, of interest to note the following remark in Ref. p|: 

'In Fig. 1.2-1, we show the two dimensional histogram of carbon energy 
versus alpha-particle energy. The contribution of the broad 1~ state at 9.6 MeV 
in 16 is the dominant feature of the histogram, while the true low energy events 
are visible in a cwn>e(i|45] group roughly along the diagonal line.' 

Certainly, in the measurements of Refs. 0, Q] which used approximately 10 fig carbon 
collector foils, the pulse height relation between 12 C and a events was strictly linear. The 
curvature may suggest that more mass per unit area than stipulated in Ref. possibly 
caused by carbon deposition during bombardment, was present in the collector foil, or that 
the data acquisition had some other problem. It is not clear what measures, if any, were 
taken to prevent carbon build-up on the foils, or if the foil thickness was otherwise monitored. 
In any case, accurate energy calibration of a curved group in the two dimensional spectrum 
appears difficult. The method of energy calibration was, however, not reported (see Sec. 
IIV A 21) . No coincidence efficiency measurement is described. 

The actual geometry, details such as beam spot size and detector distances, of the Seattle 
experiment are not known. Therefore only a simulation using very generic features can be 
performed, see Sec. IIII D[ 
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D. The Argonne measurement 



Recently a paper was published 8j in which a measurement of the /3-delayed a-spectrum 
of 16 N at Argonne National Laboratory is described. In this experiment, a fast 16 N beam 
produced by the d( 15 N,p) 16 N reaction was slowed down in a gas cell and captured onto a 
target carbon foil with a claimed beam spot of 5 mm. Within 15 s, the collected activity 
was rotated into an ionization chamber where both the a-particle and the 12 C recoil were 
detected. A similar ratio cut as in Ref. [2| has been applied, eliminating degraded target and 
detector events. Because of the particular geometry, events hitting the target frame were 
also visible in the two dimensional pulse height spectra and were eliminated. 220,000 events 
were published. The energy calibration used the 10 B(n,a) 7 Li and the 6 Li(n,a)t reactions, 
with the lowest energy calibration line at 1.472 MeV a energy. The electronic coincidence 
efficiency was checked with a pulser; however, no check of the coincidence efficiency for the 
full system has been presented. The energy resolution is quoted with 40 keV at 1.4 MeV; 
no response function of the detector is reported. 

III. GEANT4 SIMULATIONS OF /3-DELAYED a-DECAY SPECTRA OF 16 N 
A. Introduction 



For extended Monte-Carlo (MC) simulations, the GEANT4 [19| toolkit has been used. It 
provides excellent possibilities for the definition of complex geometries and particle tracking 
therein, but lacks in support for processes often relevant for nuclear physics. In particular, 
there is no process for Coulomb scattering of ions at low energies, and the energy loss 
tables are not up-to-date. The built-in multiple scattering process is only useful for high- 
energy light particles scattered on heavy ions, neglecting the energy transfered to the recoil 
particle. Therefore, a custom-made Coulomb scattering process has been developed and 
added to GEANT4 to test primarily ion propagation in ionization chambers. This process 
has been applied here. 

GEANT4 asks in each process for the mean free path, and takes the smallest of these 
values, or the distance to the next geometrical boundary, applies all "continuous" processes 
(only the energy loss in our case), and finally all "discrete" processes (only Coulomb scat- 
tering in our case). The custom-made Coulomb scattering process as well as the energy loss 
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process are described in the following. 



1. The Coulomb scattering process 

The code is based on a paper from Moller, Pospiech, and Schrieder |2(|, describing the 
implementation of a Monte Carlo (MC) simulation for a chain of single scatterings on a 
screened Thomas-Fermi potential. Using reduced energy and reduced cross-section quanti- 
ties the problem can be simplified to the integration of one scalar function (see. Ref. [21[ for 
the tabulated function) for all charges, masses, energies, and angles. In this approach, the 
total cross section is a free parameter; the larger the cross-section, the more scattering will 
be calculated, and the result will be more accurate for small scattering angles. The physical 
limit is the half-distance tq between the atoms of the material the particle is moving in, 
r = 0.5 n -1//3 , where n is the atomic density of the material. A lower limit is given by the 
screening radius a = 0.8853 as/ ' \J Z^ 3 + Z^ 3 , with as as Bohr radius. The reduced total 
cross-section J tot is defined as J to t = o /na2. For the calculations here we used values of J to t 
= 1-10. 

For each scattering, the change of direction and energy is calculated and the recoil ion is 
generated if is above a threshold (typically 1 keV). GEANT4 also tracks the recoils, so they 
obey the same processes as defined for the primary ions. 



2. The energy loss process 



The stopping power calculation in GEANT4 is based on the Ziegler parameterization [23] 
for stopping of protons in all matter, and on the effective charge model to scale the energy 
loss from protons to higher charged ions. However, different sets (tables) of coefficients can 
be selected based on different evaluations. The default set is based on the ICRU 49 report 



24| but available are also tables from Ziegler 1977, Ziegler 1985, and SRIM 2000. For the 
ICRU and Ziegler tables two variants can be selected, based on helium stopping powers and 
on proton stopping powers. For the calculations here, we used the ICRU-49 table based on 
He stoppings but we also compared the results with calculations based on other tables. The 
differences are not significant. 
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FIG. 7: (color online) GEANT4 geometry used to simulate the Mainz experiment The figure 
shows the gas cell and the outside detectors, a-particles and 12 C nuclei are multiple scattered by 
the gas, the collodion foil, and in the gold layer of the detector surface. 
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FIG. 8: (color online) GEANT4 simulations and comparison with the Mainz spectrum. The 
simulated spectrum is shown with the x-points, the simulated, shifted (63 keV) spectrum is shown 
by the cross points (+), and the experimental Mainz spectrum by the star (*) symbol. The latter 
two spectra overlap almost perfectly. The energy is in the laboratory system. 



B. The Mainz spectrum 

The exemplary description of the Mainz experiment presented in Ref. [U] allows one 
to simulate this measurement with a high degree of confidence. The geometry used in the 
GEANT4 simulation is shown in Fig. [7J The fiducial volume, i.e. the volume from where a 



particles can reach the detectors is identical to the one of Ref. 



Ill- 



The results of the simulation are shown in Fig. [HJ The simulations (using the TRIUMF 
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initial distribution) show that the Mainz spectrum would be shifted downward by about 
60 keV, if an intrinsic calibration had not been applied. Applying such a shift shows that 
the experimental and the simulated Mainz spectrum show excellent agreement. The shift 
and broadening relative to the TRIUMF spectrum confirms that any notion of a 'zero mass 
spectrum' attributed to the Mainz measurement in Ref. [3( is not justified. 



C. The Yale spectra 



1. GEANT4 Simulations of the Yale spectra 

The information available about the Yale experiments was used to simulate them. A 
16 N source was uniformly distributed through a 180/ig/cm 2 aluminum foil, at first without 
any extension in area. All angles 9 were simulated, where 9 is the initial emission angle 
of the a particle relative to the foil. (Zero degrees is perpendicular to the foil.) The 
spectrum derived from the TRIUMF R-matrix fit 2] was used in most studies for the original 
distribution, excluding thus the detector resolution in the measured TRIUMF spectrum. 
Also the final spectrum of Ref. |6J was employed once. The result of our kind of simulation is 
equivalent to a convolution of the initial spectrum and always leads to a 'broader' spectrum. 
Thus it is pointless to use as an initial spectrum one which is already 'broader' or of the 
same width, as no agreement will result. Fig. [9] shows a comparison between the original 
input spectrum, the spectrum subjected to a thick catcher as in Fig. H] and the GEANT4 
calculations for a detector opening angle of 10°. It is obvious that the GEANT4 simulated 
spectrum agrees very well with the previously simulated spectrum (Fig. H|) that includes 
only catcher thickness effects. However, at the low energy side, a tail develops caused by 
small to large angle scattering obscuring the physical low energy plateau of the 16 N spectrum 
and eventually leading to many very low energy events. Those would typically be below the 
detection threshold, or obscured by (3 rays from the abundant 16 N decays. 

The spectrum deteriorates with increasing opening angle of the detectors. This is shown 
in Fig. [10] displaying different angle ranges. As can be seen, there is a target thickness 
effect by projection, making the main peak wider as 9 increases. In addition, the Coulomb 
scattering becomes more and more significant, until at high angles, it obscures the underlying 
spectrum completely. The total counts in the spectra of Fig. [JU] scale with sin 9 until at 
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FIG. 9: (color online) Comparisons of the original spectrum (bars) used in the simulation (TRI- 
UMF), a spectrum with energy loss only applied (asterisks), and one with full scattering (crosses) 
from GEANT4 for a point source extending down to zero energy. The energy is for a-particles in 
the laboratory system. 
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FIG. 10: (color online) The GEANT4 simulated spectrum for different angular ranges of 8, as 
indicated by the labels for the Yale conditions. 

angles close to 90° the foil itself obscures detection. 

In a subdue,, St ep, we nave seated the „ of the FVanee III expert flB, 
and, as the geometry is shnilar, the Zhao expe^entS 

, including the extended catcher 
foil and the large-area silicon detector array. We have derived both single energy response 
functions for each of the energy points in the France III spectrum, as well as an integrated 
spectrum from the TRIUMF distribution. Fig. [TT] shows the single energy response function 
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FIG. 11: Response function using the Yale setup for a monoenergetic a-source at Eq,=1.71 MeV 
showing the energy loss range (plateau), a short range tail (down to about 1.2 MeV), and a long 
range tail with a many very low energy events. The energy is in the laboratory system. 
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FIG. 12: (color online) GEANT4 simulation using the Yale set-up for an a-spectrum corresponding 
to the TRIUMF distribution. 

for E Q =1.71 MeV, close to the main peak of the spectrum. As is visible, the response function 
shows the typical broadening by the catcher foil energy loss (plateau) with a short range 
tail going down to about 1.2 MeV, followed by a long range tail leading to many very low 
energy degraded events. 

In Fig. [12] a simulation is shown which uses the geometry of the Yale experiments and the 
TRIUMF a-particle energy distribution. No detector energy resolution has been included 
into the simulations which may partially explain the slightly wider than simulated France 
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FIG. 13: (color online) Spectral response by response functions taken at single energy points (see 
Fig. [TTj) employing the Yale set-up for an a-spectrum corresponding to the TRIUMF spectrum 
(bar) and the final France III spectrum (asterisks) in comparison with the France III /3-corrected 
(cross) and /^-uncorrected spectra (square). The simulation spectra have been folded with a 55 
keV Gaussian detector response 



Ill spectrum in the high energy region. We also suggest that the /3-efficiency correction in 
the high energy region of Ref. [6( is somewhat overdone (see Fig. [2]). The simulation tail is 
slightly higher than the measured one. This might be corrected by assuming a nonisotropic 
or nonhomogeneous distribution of 16 N in the catcher foil, but there is no information about 
such a possible distribution. 

As discussed above, response functions for each energy in the France III spectrum were 
calculated. From these similar spectra to that shown in Fig. [T2] can be derived for both 
the TRIUMF and the France III final spectrum. These are shown in Fig. [13] in comparison 
with the France III raw spectrum and the non-/?-efficiency corrected spectrum (see also 
Fig. [2]). As in Fig. [121 Fig- [TBI shows the simulated spectrum with the TRIUMF initial 
spectrum to be between the /3-corrected and /9-uncorrected spectra in the high energy region, 
while the France III spectral input leads to points above the corrected spectrum. However, 
the agreement with the France III spectrum at the low energy side, where /^-efficiency 
corrections are negligible, is excellent for the case of the TRIUMF spectral input, while for 
the France III final spectrum as input there is a considerable discrepancy. The simulations 
also demonstrate that the six lowest energy points in the France III raw spectrum are not 
related to the plateau in the spectrum at low energies found at TRIUMF, but to the tail in 
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the response function mainly caused by Coulomb scattering in the thick catcher foil. 

In general, the simulations show that the measurements at TRIUMF and Yale are in 
reasonable agreement (for both France III and Zhao), while the final France III spectrum is 
in disagreement with the original spectrum. 

In the subsequent section, we investigate, by applying a deconvolution method described 
in literature, the degree to which useful information can be extracted from the Yale mea- 
surements. 



2. Deconvolution of the Yale spectra 

Deconvolution of experimental information is, in general, a difficult problem. There are 
several techniques known, all of them iterative, but not convergent. In the case of the (3- 
delayed a-decay of 16 N, we are dealing with a spectrum, i.e. a finite number of discrete, 
positive data points, something we also expect as the final solution. For such a problem 



the Gold algorithm of deconvolution is most appropriate 25] 46]. The data analysis frame 



'ROOT' [26j developed at CERN provides an implementation of the Gold deconvolution 
algorithm that we have applied here, in particular by the TSpectrum::Unfolding() function. 
For the response function, necessary to be known in the deconvolution, the GEANT 4 
simulated spectra for each of the raw Yale data points, as discussed in Sec. IIII C 1[ Fig. [TTJ, 
are employed. 

As the deconvolution in the Gold algorithm is an iterative process, the number of itera- 
tions itself is a free parameter to be chosen by judgment as the process, in general, is not 
very well convergent. For the initial iterations the spectrum gets narrower and narrower 
while the low energy tail gets smaller. We find only a mild convergence for the high statis- 
tics points at the maximum of the /3-delayed a decay distribution of 16 N for an increasing 
number of iterations starting at about 20 iterations, while the low count regions at both 
the low and the high energy sides of the 16 N spectrum start to fluctuate with an increasing 
number of iterations. The width of the main peak converges to the TRIUMF width or to 
a value slightly below. In Fig. [TH deconvoluted spectra obtained by the Gold algorithm for 
increasing iteration steps are shown. 

From this deconvolution process applied here it can be concluded that, without any pre- 
conceived knowledge of the shape of the /^-delayed a spectrum of 16 N, only the presence of 
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FIG. 14: (color online) Deconvoluted spectra according to the Gold algorithm for increasing num- 
bers of iteration steps (n=l-1000) as indicated by the labels on the dashed lines. The France III 
final spectrum is shown as the boxes and the TRIUMF i?-matrix fit as straight line. The energy 
is in the laboratory system. 

the major peak with some reasonable width can be concluded, while any other information 
is obscured by the experimental procedure chosen in the Yale measurements. Poor energy 
resolution cannot be reversed. 



D. Simulations of the Seattle and TRIUMF experiments 

Using the available information, we attempted to simulate the Seattle experiment. As 
;his spectrum has used 12 C+a coincidences it was first simulated, how, despite the ratio cut 



4=71 ] . the spectral form changes with increasing target foil thickness. What happens for a 
10/ig/cm 2 C foil in comparison to a 100 /ig/cm 2 C foil is shown in Fig. [T3 The spectrum 
of the 100 /ig/cm 2 C foil has been shifted higher in energy to match that of the 10 /ig/cm 2 
C foil. In both simulations, the spectra use the TRIUMF input. It is clear that while tail 
events stay relatively well preserved, there is a short range tail component that is not easily 
removed by the ratio cut. Therefore also the method of coincidence ratio cuts calls for as 
thin foils as possible. 

In Fig. [16] a comparison between the TRIUMF based simulated spectrum based on a 
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FIG. 15: (color online) Comparison of a simulated spectrum using a 10 /ig/cm 2 C foil (x) and a 
100 fig/ cm 2 C foil (bar). The energy is in the laboratory system. 
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FIG. 16: (color online) Comparison of a simulated spectrum using a 20^g/cm 2 C foil (cross) and 
the Seattle spectrum (bar). The energy is in the laboratory system. 

20/ig/cm 2 C foil and the Seattle spectrum is shown. A slight energy adjustment, corre- 
sponding to the energy loss in the foil has been done. Both the high and the low energy side 
of the Seattle spectrum are wider than in this comparison, see also the discussion in Sec. 
II V A 21 regarding energy calibration. It may be noted that approximate agreement on the low 
energy side can be obtained with a foil thickness of approximately 50 /zg/cm 2 . However, no 
possible foil thickness in the simulation can reproduce the high energy side of the spectrum 
from the TRIUMF input. 

As seen for the simulations to the Seattle data a 10 /xg/cm 2 C foil introduces little 
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FIG. 17: (color online) Simulated two dimensional energy versus energy spectrum for a gas chamber 
with 150 Torr as described for the Argonne experiment. The asymmetry in regard to the 45° axis 
originates from the asymmetric collector foil mounting in reference to the foil supporting frame. 
The energies are in the laboratory frame. 

disturbance to the /3-delayed a spectrum of 16 N. We have used the TRIUMF geometry and 
input to simulate the TRIUMF data as well. The agreement is good after folding with the 
approximately 30 keV energy resolution of the detectors. 



E. Simulations of the Argonne experiment 

Following the description of the experiment 8|, the spectra have been simulated in 
GEANT4 based on the TRIUMF spectrum as previously. No 16 N isotopes in gaseous form 
and no noise, as reported in Ref. 8j, have been included. The two-dimensional energy versus 
energy spectrum from the simulations is shown in Fig. [T71 Particles were tracked through 
the P10 gas (90% argon, 10% CH 4 ) to complete stop either in the gas or a boundary. The 
reported partial energy collection for high energy events for ion chamber pressures of 150 
Torr [8] are confirmed. 

In Fig. [18] the comparison between the simulated GEANT4 spectrum and the spectrum 
of the Argonne work is shown. In the simulated spectrum a condition (cut) for maximum 
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FIG. 18: (color online) Comparison between the simulated (bar) and the Argonne spectrum (x). 
A pressure of 195 Torr is used. The energy is the a-energy in the laboratory in keV. 

carbon energy (channel) as well as a ratio condition were applied. As the numerical value 
for the latter is not reported in Ref. [sj], we applied one similar to Ref. Q]. The figure shows 
reasonable agreement on the low energy side of the main peak, with a slight excess of counts 
in the simulated spectrum around E a =1500 keV and E a = 1 MeV, see also the discussion 
in Sec. IIVA4I The high energy side of the spectrum is definitely wider, but no detector 
resolution effects have been included in the simulation. The simulation demonstrates also 
that measurements using ionization chambers have potential problems with the response 
tail from the target, as this limit seems to be reached in the Argonne measurement. In 
conclusion: the Argonne measurement is consistent with the TRIUMF one at least for the 
region of E a =l.l to 1.3 MeV. For further discussion, see Sec. IIV A 41 



IV. COMPARISON AND R-MATRIX FITS TO DIFFERENT DATA 



A. Fits and comparison of spectra 

Independent of their credibility all of the /3-delayed a spectra of 16 N can be subjected 
to R matrix fits. In particular, this can show how much the apparent differences between 
the spectra influences the deduced S-factor S^i(300). Only fits with an interaction radius 
of a=5.5 fm and with the j3 feeding factors A23=0 and ^33=0 for the /-wave (as in Ref. 
also see Ref. [2| for i?-matrix notation) were performed, as it is not the goal of this paper 
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FIG. .19: (color online) Comparison of (he Una] j-dda.yed n-spertra of Ref. | [j], [f| (bars) and 
Ref. [^J (crosses) with both spectra normalized in height to the one of Ref. Q]. The energy axis is 
the center-of-mass-system. Also shown are i?-matrix fits to these spectra. 



to produce a final combined result with a variation of all parameters. However, it must be 
stressed again that the results derived from meaningless spectra are meaningless. 

Comparisons of the final experimental spectra are as long meaningless as all of the mea- 
surements have different energy resolutions and different 16 N source thicknesses. In principle, 
only unconvoluted .R-matrix fits, with convolutions applied in the fittings, should be com- 
pared, as clearly a spectrum with a detector energy resolution of 10 keV will be different 
from one with a resolution of of 100 keV, even, if the detector responses would be simply 
Gaussian over many orders of magnitude. 



1. Comparison between the France III and the Zhao spectra 

In Refs. no comparison between the 16 N spectrum presented there and the one 

of Ref. is given, even though differences with the spectrum of Ref. [^j are pointed out 
with considerable detail. Both spectra , together with R-matrix fits to them, are shown 
in Fig. [191 Several observations may be made: the two spectra do not agree, not only on 
the low energy side of the main peak, but also on the high energy side even though they 
both originate from the same initial distribution. The quality of the i?-matrix fits for both 

n 

spectra is more than excellent, as already noted above for the case of Ref. [3|]. Because 
both spectra have gone through extensive, non-justified, data reduction, it is also unclear, 
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(t. 20: (color online) Comparison of the (inal . /-delayed a-spectra of Ref. ?J (crosses) and Ref. 
(bars) with the Seattle spectrum normalized to that of Ref. at the peak. Energies are in the 
center-of-mass system. Also shown are ii-matrix fits to these spectra. 



how the detector resolutions should be treated. E.g. for the Mainz spectrum [ll], a detector 
resolution of 12-18 keV (FWHM) at 1.5 MeV is quoted, with the full system response leading 
to about 40 keV at 1.2 MeV a energy, while for Ref. 15|] 50 keV is given at an unspecified 
energy. However, as discussed above, the spectrum of Ref. 3| has been normalized to that 
of the Mainz group. Indeed, a low energy resolution (20 keV or less) fits the spectrum of 
Ref. [3j] best, particularly around the maximum. On the other hand, the broad spectrum 
of Refs. 14 a ]& is well fitted with the quoted experimental detector resolution of 55 keV. 

n n r 

The particular fits shown in Fig. [T9l correspond to S^i(300)=79.6 keV b for Refs. [4j, |5j, 6] 
for a bad resolution fit and S^i(300) = 101 keV b for Ref. |3|] for a low resolution fit. Ref. 3] 



suggests S^i(300)=95 keV b, see also Sec. Ill B 21 No value of S^i(300) has been suggested 
for the France III 6|] measurement. 



2. Comparison between the Seattle and the TRIUMF spectra 

In Fig. [20] a comparison between the Seattle and the TRIUMF spectra are shown . 
Obviously, the two spectra disagree, as the Seattle spectrum is wider both in the low and the 
high energy regions of the TRIUMF spectrum. However, as the method of energy calibration 
for the Seattle spectrum is unknown to us, it may be noted that a simple linear recalibration 
would give a nearly perfect match between the two spectra (E new = 0.95*^^ + 0.118 MeV). 
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Note the following quotation from Ref. 181 ]: 



'Our results are consistent with the previous measurements at TRIUMF Q] 
and Yale QQ'. 

This suggests that the Seattle group does not consider the differences between those spectra 
as important. In fact, the difference between these spectra and the one of France III is the 
most pronounced. 

The detector resolution in the Seattle experiment is not published. In principle, a larger 
energy resolution can achieve agreement between the TRIUMF and the Seattle spectra. We 
have varied that resolution and find an optimum value at about 30 keV with Sei (300) =97.1 
keV b. However, fitted simultaneously with the radiative El ground state data, a smaller 
resolution (20 keV) is preferred with S^i(300)=95.7 keV b. For larger resolutions than 30 
keV, S^i(300) declines with a decrease in fit quality. No S£i(300) has been suggested by 
the Seattle group. 

3. Comparison of the Mainz and TRIUMF spectra 

A comparison between the Mainz single count spectrum and the TRIUMF coincidence 



f. Q 



jeen claimed the two both spectra 



that the TRIUMF singles spec- 



spectrum is shown in Fig. 15 of Ref. [2|]. It has never 
agree. However, it was found and pointed out in Ref. 
trum, before application of the coincidence cuts, agrees well with the Mainz spectrum. The 
difference with the coincidence spectrum is largely from the coincidence ratio cuts applied 
to remove degraded a events found in the coincidence spectrum. In the simulations of Sec. 
IIII Bl it has indeed be shown that the Mainz (single detector) spectrum is fully consistent 
with the TRIUMF one. 

4- Comparison between the TRIUMF and Argonne spectrum 



In Ref. [8j a comparison between the low energy part of the TRIUMF spectrum and 
the Argonne spectrum as well as with the Z. Zhao spectrum jjjj], is presented. Better agree- 
ment of the Argonne spectrum with the Mainz and Z. Zhao spectrum is claimed than with 
the TRIUMF spectrum. The comparison is in principle invalid as the TRIUMF combined 
detector and f3-v recoil resolution is 30±5 keV, while a resolution of 40 keV (no error) at 
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FIG. 21: (color online) Comparison of the /^-delayed a-spectra of 16 N Ref. |2| (bars) and Ref. 
(cross). 



1472 keV is claimed in Ref. 
43 keV 



8j ; combined with the (3-v recoil resolution that averages to 
491 ] . Fig. |2~T1 shows a comparison with the TRIUMF spectrum 501] . 
We make the following observations: As described, the region of the Argonne spectrum 
from about E Q =1.05 MeV to about 1.3 MeV is above the TRIUMF \l\ spectrum. However, 
there are other differences: (i) the high energy points of the Argonne spectrum are above the 
TRIUMF spectrum starting at about E Q =1900 keV. We attribute this to the poorer detector 
resolution in the Argonne experiment (see arguments above for non-compatibility of spectra) 
and also 18 N background effects, (ii) Also in the region near E a =1.5 MeV, the Argonne data 
are above the TRIUMF data. In that region, the continuation of the response function that 
cannot be removed by the ratio cut in the TRIUMF spectrum, has been subtracted, see Fig. 
10 of Ref. 2j. No removal of a similar response function in the Argonne spectrum, nor of 
their additional background from 16 N in the gas phase have been reported. We also note 
that points in the Argonne spectrum from about E Q =850 keV to 1050 keV are systematically 
higher than the TRIUMF points. This, however, can be resolved by different normalizations 
as the main peaks do not need to match perfectly. 



n Fig. [22] we show 
281 ] . No particular 



Most interesting are the two points at about E Q =1074 and 1094 keV. 
the previous Figure [21] with the scaled /3-delayed 18 N spectrum of Ref. 
adjustment in the energy scale of either spectrum has been done. The agreement with the 
main peak of the 18 N spectrum for these two points is remarkable. The lower energy point 
is about 3.6a above its expected value and the higher energy point 2.6a. The statistical 
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FIG. 22: (color online) Comparison of the /3-delayed a-spectra of 16 N Ref. (x), Ref. (bar) 
with the scaled 18 N spectrum of Ref. [28] (asterisks). 



deviation has been determined from R-matrix fits excluding this region, see Sec. IIVBI The 
random chance for such a peak at the energy of the 18 N main peak is about 3xl0 -6 . This 



strongly indicates a background of N in the beam. While it is argued in Ref 



that 



the combination of 15 N+d cannot produce any 17 ' 18 N, fusion evaporation reactions of the 
15 N beam with carbon or other light elements at windows or collimators can produce both 
isotopes. 

While 18 N is a relatively trivial contamination (see Sec. IIV C 21) showing mainly up as one 
narrow peak and being removable by decay time cuts, 17 N is not. The /3-delayed a spectrum 



of 17 N 29| has no narrow peak, and the half life (4.2 s) is near that of 16 N (7.1 s). Removing 
both the likely 18 N background and an 17 N background of roughly the same intensity leads 
to the spectrum shown in Fig. [231 The spectrum derived by such a background subtraction 



agrees well both on the plateau and in the region above up to E a =1300 keV with the 
TRIUMF spectrum. The low yield point at about E a =1054 keV is statistically consistent 
with the TRIUMF spectrum. In its calculation, the method of the subtraction, and details 
like matching energy calibrations and different detector resolutions that are involved pose 
rather subtle problems that are not worth improving at this stage. 
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FIG. 23: (color online) Comparison of the /3-delayed a-spectra of 16 N of Ref. |8| (asterisk) with the 
one with 17,18 N background removed (cross), as described in the text. In addition, the TRIUMF 
spectrum is shown. 



B. i?-matrix fits to the Argonne spectrum 



The Argonne spectrum has been fitted at different interaction radii a. In general, we find 
the S-factors S^i(300) to be smaller than given in Ref. 8]. However, as our spectrum is 
read out from a paper, we encounter large deviations from the fit at the peak of the 16 N 
distribution due to the difficulties of reading out the spectrum there. While it is possible to 
adjust the points to lower the least squares sum considerably, this may lead to biased fits. In 
Fig. [IM]we show a fit done with a=5.5 fm. For the resolution, we chose 57 keV in the center- 
of-mass, very close to the optimum (minimum \ 2 ) setting. For this particular fit Ssi(300) is 
found to be 68 keV b. As it is not described in Ref. [8J what constitutes a statistical error, 
we cannot derive those. In addition, it is not mentioned what fit parameters have been used, 
and what degrees of freedom are included in the partial least squares parameters. 

The Argonne group has put great emphasis on the region between Eq,=1.0 to 1.3 MeV 
showing a deviation from the TRIUMF data. To see the relevance of this region, the data 
points there have been taken out of the fit. Fig. [25] shows the fit in the region of interest 
with the data used and not used in the fit. Obviously, when those data are excluded, the fit 
runs close to the TRIUMF data. For the same exercise with the TRIUMF data the fit does 
not change significantly. 

The deviation of those points identified as likely 18 N events was taken from the fit that 
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FIG. 24: (color online) it!-matrix fits to the Argonne spectrum with a=5.5 fm. Shown are the 
spectrum (bar), the resolution convoluted fit (long dash), the unconvoluted total fit (short dash), 
the p-wave component (dotted) and the f-wave component (dash-dotted). The energy is in the 
center of mass system. 
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FIG. 25: (color online) i?-matrix fits to the Argonne spectrum with a=5.5 fm for the energy region 
displayed. Shown are the spectrum (bar), the resolution convoluted fit (long dash), and the fit, 
with those data excluded (short dash). The energy is the center of mass energy. 



had those events removed, see Sec. IIV A 41 
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FIG. 26: (color online) Comparison of the /3-delayed a-spectra of 16 N Ref. (bar) and Ref. 
(as energy-calibrated in Ref. [6J) (asterisks) with the spectrum where previously subtracted events 
from the /3-delayed a-decay of 17 N have been added back (cross). [Only the relevant low energy 
region is shown.] The energy is the center-of-mass energy. 



C. Background subtraction in the TRIUMF spectrum 



It has been suggested 



, |7| 5l| that the final TRIUMF result is questionable beyond the 
errors quoted, because a background of events from the /3-delayed a-decay of both 17 N 29] 
and 18 N [28] has been subtracted to obtain the final spectrum. While the total number of 
these subtracted events is only 2130 compared to 1,026,500 total events after subtraction, 
it has been implied that those subtractions, if reversed, would (i) reconcile the differences 
between different measurements, and (ii) cause changes larger than the assigned errors given 
to the final S-factor S^i(300) quoted in Ref. [2]. We show below that these arguments are 
without merit, and could have been checked by the authors of Refs. themselves, as 

both the total number of events subtracted as well as the shape of these backgrounds are 
given in Ref. [2], particularly in Fig. 10. 



1. The subtraction of the (3-delayed a-spectrum of N 

In Fig. [26] we compare the final TRIUMF 16 N spectrum, the TRIUMF 16 N spectrum with 
all 17 N events added back, and the Mainz spectrum (as energy calibrated by the Yale group 



6j). Clearly, the effect of the subtraction of the 17 N spectrum from the TRIUMF spectrum is 
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minor, and does not produce any reasonable identity with the Mainz single event spectrum. 
The S-factor at 300 keV is found to be S B i(300)=82.2 keV b with the 17 N events added back 
to be compared with S#i(300)=79±21 keV b with the 17 N events subtracted 



in Ref. 



21 Note that 



2| an error of ±5 keV b is given for the systematic uncertainty in the R-matrix fits 
resulting from the 17 N event subtraction. It is therefore concluded that a claim of significant 
uncertainties resulting from the subtraction of 17 N events is irrelevant. 

2. The subtraction of the (i-delayed a-spectrum of 18 N 

The /^-delayed a emitter 18 N has a half life of 0.63 s. As reported in Ref. 2|, the 
implantation time for the radioactive 16 N beam was 3 s with 0.25 s moving time between 
three subsequent decay stations. While there is small but clear evidence of the low energy 
1.081 MeV peak from 18 N in the spectra from the first detector station, this peak is no 
longer visible in the subsequent two stations as expected from the short 18 N half life. The 
four later spectra were then used to remove the 18 N contamination from the 16 N spectra. 
This information has been given previously in Ref. Q]. 

The low energy 18 N peak allows for an easy scaling of the 18 N subtraction. In Fig. [27] 
we compare the TRIUMF 16 N spectrum, the TRIUMF 16 N spectrum with all 18 N events 
added back and the Mainz spectrum, as calibrated by the Yale group [6J]. Again, it is 
clear that agreement with the Mainz spectrum is not significantly improved by adding 18 N 
events back. The fit to the added-back spectrum has a relatively poor quality as it cannot 
accommodate the narrow 18 N peak. However, leaving the 18 N events in the spectrum leads 
to an Siji(300)=74.2 keV b, still well within the quoted error of Ref. [2j, i.e. 79±21 keV 
b. Note, however, that the 18 N a peak in the 16 N spectrum provides an excellent intrinsic 
a-particle energy calibration with a precision of about 1 keV in this energy region. 

V. CONCLUSION 

In Table I we give an overview over the 16 N decay experiments discussed here. Discussing 
the rows of Tab. I: Aside from the TRIUMF experiment, all others use the 15 N(d,p) 16 N 
reaction for 16 N production. However, only the Mainz experiment 



lfl, the TRIUMF and the 



Argonne one thoroughly separated the region of 16 N production from the detection region 
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TRIUMF-backadded 



E(MeV) 



FIG. 27: (color online) Comparison of the /3-delayed a-spectra of N Ref. [2| (bars) and Ref. 
(asterisks), (energy calibrated in Ref. [6]), with the TRIUMF spectrum to which all events from 
the /3-delayed a-decay of 18 N have been added back (cross). [Only the relevant low energy region 
is shown.] The energy is in the center-of-mass system. 



TABLE I: Overview of the /3-delayed a decay experiments discussed in the article. 



Experiment 


Mainz W 


TRIUMF [2J 


Yale 1 [3J 


Seattle [7] 


Yale 2 [6] 


Argonne £8j 


16 N production 


15 N(d,p) 


Isotope Separator 


15 N(d,p) 


15 N(d,p) 


15 N(d,p) 


15 N(d,p) 


1<5 N implantation speed 


low 


low 


high 


high 


high 


high 


mass/ composition 


30 Mg/cm 2 C 6 H 7 N 2 .50io 


10 fig/cm 2 C 


180 (ig/cm 2 Al 


10/20 ng/cm 2 C 


ISO fig/ cm 2 Al 


17 Aig/cm 2 C 


of 16 N catcher 


+ 0-1.5 cm 6-8 Torr 15 N 2 gas 
cqv. 0-17 /ig/cm 2 N 












detectors 


<35 M m Si 


10.4-15.8 pm Si 


50 (im Si 


15-20 pm Si 


50 /im Si 


ion chamber 


background (measured) 


/3-tail 


1S N, 17 N 


none 


unknown 


none 


maybe 


degraded event suppression 


none 


a - 12 C 


/3-a 


Q - 12 C 




12 P 
a- C 


efficiency corrections 


none 


none 


P 


u nknown 


& 


none 


energy calibration 


10 B(n,a) 


18 N, 20 Na 


10 B(n,a) 


unknown 


10 B(n,a) 


10 B(n,a) 


dcconvolution applied 


none 


none 


division (sec text) 


unknown 


division (sec text) 


none 



While care may have been taken to avoid 16 N hitting e.g. the catcher foil frames in the Yale, 
Seattle and Argonne experiments, a deep implantation of 16 N ions has not been excluded at 
the level of 1/10,000 leading to potentially degraded a events. 

While the Mainz spectrum has been dubbed a 'zero thickness spectrum' in Ref. jjj, it 
is certainly not the one that uses the least amount of material to hold the 16 N. As the 
GEANT4 calculations show, the combination of an extended volume source and the foils in 
front of the detectors leads to some broadening. In contrast, the 10 /ig/cm 2 C foils in the 
TRIUMF experiment slowly decreased in thickness in the course of the implantations, until 
a hole was produced by the stable beam at mass A=30, and the four active 16 N collector 
foils had to be replaced. Spectra collected at different degrees of foil degradation showed no 
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difference. From Ref. 11] it is indeed known that for the Mainz spectrum the accumulated 
mass traversed by the a-particles leads to an increase from about 15 keV detector resolution 
to 40 keV total resolution. The Mainz spectrum also allowed a relatively large angular range 
of collection, leading to a slightly increased range in matter. As far as the Seattle spectrum 
is concerned, it is not clear whether the spectrum that has been distributed is from a 10 or 
20 /ig/cm 2 foil It also can not be excluded that additional carbon build up was involved. Of 
course, the catcher foils used by the Yale group are much the thickest and, as shown above, 
this leads to heavily distorted spectra dominated for low energies by scattering. 

As far as detectors are concerned, all groups but one use silicon based detectors, with the 
thinnest detectors, i.e. those with the least (3 energy losses and thus f3 response used by the 
TRIUMF group. The thickness of the gas layer used to stop the a-particles in the Argonne 
experiment is comparable with the TRIUMF detectors. 

Both, in the Mainz and the TRIUMF spectra, backgrounds have been subtracted. This 
issue is extensively discussed above for the TRIUMF spectrum (Sec. IIV CI) . In fact, for the 
Mainz spectrum, the (3 induced background from the 16 N decay is quite significant at low 
energies. 

As the primary goal of the Mainz experiment was not to determine the shape of the 16 N 
spectrum exactly, but to search for the narrow line of the parity- forbidden a-decay of 16 O 
at 1.283 MeV of a energy, no measures were taken to check for degraded events. However, 
in all other experiments, some measures for increased discrimination against degraded or 
otherwise unwanted events were taken. In the case of the Yale experiments (3-a coincidences 
were employed. Those, in principle, eliminate events from the (3 tail of the a-detector and 
from events, where a particles hitting the detector do not deposit the full charge. However, 
no means exists in this method to eliminate events that suffer from some degradation of 
a-particle energies in the thick aluminum foil. This is consistent with the observations and 
simulations about these measurements discussed above. As the detection of very low energy 
/3-particles is required for high energy a-particles, extensive (3 efficiency corrections were 
necessary, in addition, to derive the two Yale spectra. 

In the Argonne, Seattle and TRIUMF experiments a- 12 C coincidences with subsequent 
pulse height ratio cuts were used. That, in principle, removes /3-tail events, a good fraction 
of degraded detector events and events scattered in the collection foil. As our simulations 
show, a low level of degraded response events will remain depending on catcher foil thickness. 
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While we do not know the method of energy calibration for the Seattle experiment, the 
Argonne, Yale and Mainz experiments employed the 10 B(n,a) 7 Li reaction for their energy 
calibration. The reaction results in two a-particle groups, at E a =1472 and 1777 keV, i.e. 
close together and close to the 16 N main peak. One also does obtain two lower energy 7 Li 
groups. However, those are less suitable for calibration, as different pulse height defects or 
other differences for the differences in the detection of the 7 Li ions are difficult to take into 
account. Therefore, no low energy calibration is available in these experiments. It should 
be noted that the method of energy calibration is rather uncritical for the Yale experiments 
as they are shifted upward by about 200 keV and use an invalid deconvolution method. 
In addition, Argonne used the 6 Li(a,n)t reaction for a high-energy calibration point. The 
Argonne measurement claims a smaller systematic error from their energy calibration (5%) 
than the TRIUMF one (10%) Q. 

While not having employed it as a calibration in the TRIUMF experiment, we have 
observed that the actual position of the a peaks from 10 B(n,a) 7 Li was not consistent with 
the calibration obtained from 18 N and 20 Na. Indeed the calibration with 10 B(n,a) 7 Li was 
found to be dependent on the direction of the incident neutron flux, suggesting that only 
incomplete neutron thermalization was achieved. Complete thermalization may be hard to 
prove without an accurate measurement of the neutron spectrum at the position of the 10 B 
target. 

It is our opinion that the best energy calibration comes from the narrow decay lines 
of well known /3-delayed a-emitters that can be used without modifying the experimental 
set-up. 

Besides the Yale measurements, no other measurement relies on deconvolutions. It has 
been demonstrated above that both methods applied are incorrect. It also has been shown 
that both original measurements are consistent with the TRIUMF data when compared with 
simulations and, for the case of France III |6j, inconsistent with their own 'deconvoluted' 
spectrum. 

The TRIUMF spectrum, among all measurements, shows the narrowest main peak. It 
has been shown above that the energy calibration used for the derivation of this spectrum 
is the best one available, and that the shape of the TRIUMF spectrum is little influenced 
by background subtractions. In principle, besides resolution broadening, also an energy de- 
pendent efficiency can lead to a broadening, and more likely narrowing (efficiency dropping 
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off at low energies) of the spectrum. However, it is shown experimentally in Ref. |2j that 
the coincidence detection efficiency stays constant over the energy range covered in the mea- 
surement. None of the other measurements considers this point. As far as broadening by 
target and detector effects is concerned to the authors of this article there is no reasonable 
physical response function known that would transform the wider spectra of the other mea- 
surements into a narrower one, as measured at TRIUMF. But there are many well known 
transformations that can produce a wider spectrum from a narrow one. 



VI. OUTLOOK 



We have presented fits to several spectra in this article. However, either the derivation 
of those spectra is in doubt, even if one does not take systematic errors into account, as for 
the Yale spectra, or the spectra are not completely published, as for the Argonne spectrum. 
Therefore any derivation of a common S-factor Sei(300) would lead to questionable results. 
We therefore reiterate the value of S,ei(300) derived in Ref. [2j, i.e. S E i(300)=79±21 keV 
b includin g sy stematic errors. A new evaluation of the entire reaction rate of 12 C(«,7) 16 



is in work 



31] taking new information from other measurements into account, particularly 



phaseshift and radiative capture data. 

While other partial cross sections of 12 C(a,7) 16 have relative uncertainties exceeding 
those of the El ground state transition, in the long ton, it ce rt ain,y won,d be desirable to 
improve the value derived in Ref. |2j to a relative error of less than 10%. While statistics 
plays some role, experimental limitations and systematic errors are the real limits. 

Some numbers from other measurements are part of the fit to the /3-delayed a-spectrum 
of 16 N. These are the energy of the subthreshold 1~ state, its radiative width, and the f3- 
decay branching ratio into the E 3 .=9.6 MeV 1~ state in 16 0. While the former numbers 
appear reasonably well known, there is some uncertainty about the latter branching ratio. 



The compilation of Ref. 



321 ] quotes 1.20x10 5 ±0.05 for this branching ratio citing Ref. 



as the origin of this number. However, Ref. [33j is a theoretical work quoting Ref. 
34] . another theoretical work. Ref. 34| quotes Ref. (an early publication of the 

Mainz group) which gives a branching ratio of 1.19x 10~ 5 ±0.10. Both theoretical works 
do not quote an experimental error, so the origin of the error in Ref. 32jis unknown to 
us. However, Ref. p. 327, gives a revised value of 1.13x 10~ 5 ±0.08 for the branching 
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ratio into the E x =9.6 MeV state of 16 0, a not insignificant deviation. Therefore a revised 
measurement with improved errors of this branching ratio is desirable. 

Among the errors inherent to the measurement is the energy calibration of the spectrum. 
Ref. Q quotes a systematic error of 10 keV b for S#i(300), while Ref. 30] quotes only 5 
keV b with a more indirect measurement. This error can certainly be improved. 

Part of the uncertainties in the fits originates from the presence of both a p- as well as 
an /-wave in the /3-delayed a-spectrum of 16 N. Higher statistics data at very low energies 
(E a < 600 keV) would be very desirable in this context as they originate only from p-wave 
decay, while due to penetrability reasons the / wave drops off faster and at higher energies. 
However, this a-energy region is difficult to access, first because of possible f3~ back ground 
in the detectors that eventually may be removed, but more important second because of 
a possible response of high energy a events into this region. Our simulations have shown 
(Sees. HHP] IIIIEj) that indeed the most advanced experiments using carbon foils of about 
10 /zg/cm 2 and coincidence techniques are close to this response limit. Foils with much less 
mass seem to be impractical. As an alternative 16 N could be trapped in electromagnetic or 
opto-magnetic traps with no substrates for the 16 N involved. Because of trapping efficiencies 
and likely difficult geometries such an experiment will require high yields of 16 N. In this case 
careful simulations of the next order of response like those originating from the detectors 
are necessary. Another possibility, in principle, is to measure the (3-a correlations in the 
decay of polarized 16 N that would differentiate the data between the p- and the /-wave in 
the region of about and above E a =l MeV, see Ref. 16] . 
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and 'convolution' with a response function G(E, E — E') is made. 

We label the experiments subsequently by the laboratory where they took place. 

There is an additional error in Table A.l[6] as the points at 1468 and 1482 keV have the same 

number of counts. The correct number of counts of the point at 1468 keV should be 3785. 

[...] notes here and in other quotations an omission of text for clarification. 

Of course, if the final result is known, it is possible to construct a response function G(E, E' — 

E) = tlj^5{E' — E) with fi(E) the initially measured spectrum and ff(E) the 'unfolded' 

spectrum, see Sec. Ill B 21 Such a function, however, is not independent of the initial (a) 

spectrum. 

Highlighted here. 

It may be noted that any attempt to find an R-matrix solution based on three states for the 



/-wave presented in Ref. 



3|] leads to no acceptable result. 



The Mainz measurement [l]] certainly employed a larger catcher thickness than the TRIUMP^] 
one, but none of the experiment had zero target thickness, see Sec. IIIIB1 
There is no description how the errors were propagated, but it appears, they also went from 
original to final by the same division. 

A response function is defined by the energy dependent function (spectrum) derived from the 
propagation of events from a monoenergetic (a) source through the experimental apparatus. 
Emphasized here. 

In their introduction the authors of Ref. 



251 ] make the following statement about deconvo- 
lution: 'From the numerical point of view, deconvolution belongs to one of the most critical 
problems. It is a so-called ill-posed problem, which means that many different functions solve 
the convolution equation within the allowed error bounds. The estimates of the solution are 
extremely sensitive to errors in the input data. Very frequently the noise present in the input 
data causes enormous oscillations in the results after deconvolution.' 
[47] We call a 'ratio cut' a condition on the spectra where only a limited range of ratios of apparent 
a and 12 C events is considered to be valid. If there were no resolution and pulse height issues, 
the energetic ratio between a- and 12 C events would be exactly 3. For a description and 
application of such a condition, see Ref. [2]. 
[48] Citations as quoted here. 

[49] In the i?-matrix fits of Ref. Q] it is not mentioned, if the fits do include a convolution by the 
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detector resolution. (3-v recoil effects are not mentioned in Ref. 
publication of 16 N spectra, except the THIUMF one fl. 



and for that matter in any 



[50] As the Argonne group has not agreed to make the digitized spectrum publicly available [2^], 
we have read it out from the publication. Minor inconsistencies particularly for the high count 
points, are therefore possible. Errors are taken as the square root of the counts, as apparently 
has been done for the low count points, where errors can be read out. 
[51] Quoting from Ref. 6[:'The disagreement could quite possibly be due to over subtraction of 18 N 
contamination in the TRIUMF data. [...] The disagreement in the region of the interference 
minimum around 1.4 MeV is sufficient to change the f-wave component and leads to imprecise 
determination of the p-wave contribution to the 16 N spectrum.' Obviously an imprecision 
beyond the errors presented by the TRIUMF experimental group is meant. 
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